Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China          下载免费PDF全文

S.‐J. Wang  L. Wang  M. Brown  P. M. Piccoli  T. E. Johnson  P. Feng  H. Deng  K. Kitajima  Y. Huang 《Journal of Metamorphic Geology》2017,35(6):601-629

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.  相似文献   

Glacial flood pulse effects on benthic fauna in equatorial high‐Andean streams          下载免费PDF全文

Sophie Cauvy‐Fraunié  Patricio Andino  Rodrigo Espinosa  Roger Calvez  Fabien Anthelme  Dean Jacobsen  Olivier Dangles 《水文研究》2014,28(6):3008-3017

Equatorial glacier‐fed streams present unique hydraulic patterns when compared to glacier‐fed observed in temperate regions as the main variability in discharge occurs on a daily basis. To assess how benthic fauna respond to these specific hydraulic conditions, we investigated the relationships between flow regime, hydraulic conditions (boundary Reynolds number, Re*), and macroinvertebrate communities (taxon richness and abundance) in a tropical glacier‐fed stream located in the high Ecuadorian Andes (> 4000 m). Both physical and biotic variables were measured under four discharge conditions (base‐flow and glacial flood pulses of various intensities), at 30 random points, in two sites whose hydraulic conditions were representative to those found in other streams of the study catchment. While daily glacial flood pulses significantly increased hydraulic stress in the benthic habitats (appearance of Re* > 2000), low stress areas still persisted even during extreme flood events (Re* < 500). In contrast to previous research in temperate glacier‐fed streams, taxon richness and abundance were not significantly affected by changes in hydraulic conditions induced by daily glacial flood pulses. However, we found that a few rare taxa, in particular rare ones, preferentially occurred in highly stressed hydraulic habitats. Monte‐Carlo simulations of benthic communities under glacial flood reduction scenarios predicted that taxon richness would be significantly reduced by the loss of high hydraulic stress habitats following glacier shrinking. This pioneer study on the relationship between hydraulic conditions and benthic diversity in an equatorial glacial stream evidenced unknown effects of climate change on singular yet endangered aquatic systems. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data          下载免费PDF全文

Antonio‐Juan Collados‐Lara  Eulogio Pardo‐Igúzquiza  David Pulido‐Velazquez 《水文研究》2017,31(10):1966-1982

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.  相似文献   

Petrological evidence for shock‐induced high‐P metamorphism in a gabbro          下载免费PDF全文

J.‐J. Yang  H.‐R. Zhang  A.‐P. Chen  M.‐X. Huang 《Journal of Metamorphic Geology》2017,35(2):121-140

Microlites (minute spherulitic, dendritic, skeletal, acicular and poikilitic crystals) diagnostic of crystallization in quenched melt or glass in fault rocks have been used to infer fossil earthquakes. High‐P microlites and crystallites are described here in a variably eclogitized gabbro, the wallrock to the coesite‐bearing eclogite breccia at Yangkou in the Chinese Su‐Lu high‐P metamorphic belt. The studied hand specimens are free of discernible shear deformation, although microfractures are not uncommon under the microscope. In the least eclogitized gabbro, the metagabbro, stellate growths of high‐P minerals on the relict igneous minerals are common. Dendritic garnet crystals (<1?5 μm) grew around rutile and/or phengite replacing ilmenite and biotite, respectively. Skeletal garnet also rims broken flakes of igneous biotite and mechanically twinned augite. Radial intergrowths of omphacite and quartz developed around relict igneous orthopyroxene and are rimmed by skeletal or poikilitic garnet where a Ti‐bearing mineral relict is present. Acicular epidote, kyanite and phengite crystallites are randomly distributed in a matrix of Na‐rich plagioclase, forming the pseudomorphs after igneous plagioclase. In the more eclogitized gabbro, the coronitic eclogite located closer to the eclogite breccia, all the igneous minerals broke down into high‐P assemblages. Thick coronas of poikilitic garnet grew between the pseudomorphs after igneous plagioclase and ferromagnesian minerals. The igneous plagioclase is replaced by omphacite crystallites, with minor amounts of phengite and kyanite. Thermodynamic modelling of the plagioclase pseudomorphs shows an increase in P–T in the wallrock from the metagabbro to the coronitic eclogite, and the P–T variation is unrelated to H₂O content. The fluid‐poor pressure overstepping scenario is unsupported both by phase diagram modelling and by whole‐rock chemical data, which show that the various types of eclogitized gabbro are all fairly dry. A large pressure difference of >2 GPa between the metagabbro and the coesite‐bearing eclogites ~20 m apart cannot be explained by the subduction hypothesis because this would require a depth difference of >60 km. The microlites and crystallites are evidence for dynamic crystallization due to rapid cooling because constitutional supercooling was unlikely for the plagioclase pseudomorphs. The lack of annealing of the broken biotite and augite overgrown by strain free skeletal garnet is consistent with a transient high‐P–T event at a low ambient temperature (<300 °C), probably in the crust. Therefore, the eclogitization of the wallrock to the eclogite breccia was also coseismic, as proposed earlier for the eclogite facies fault rocks. The outcrop‐scale P–T variation and the transient nature of the high‐P–T event are inconsistent with the other existing tectonic models for high‐P metamorphism. The fact that the less refractory but denser biotite is largely preserved while the more refractory but less dense plagioclase broke down completely into high‐P microlite assemblages in the metagabbro indicates a significant rise in pressure rather than temperature. Given that the metamorphic temperatures are far below the melting temperatures of most of the gabbroic minerals under fluid‐absent conditions, stress‐induced amorphization appears to be the more likely mechanism of the coseismic high‐P metamorphism.  相似文献   

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation          下载免费PDF全文

Chongbin Zhao  Thomas Poulet  Klaus Regenauer‐Lieb 《国际地质力学数值与分析法杂志》2015,39(6):571-593

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions          下载免费PDF全文

Xiong Zhang  Jean‐Louis Briaud 《国际地质力学数值与分析法杂志》2015,39(13):1369-1409

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Soil property variation mapping through data mining of soil category maps          下载免费PDF全文

Fei Du  A‐Xing Zhu  Lawrence Band  Jing Liu 《水文研究》2015,29(11):2491-2503

Spatial information on soil properties is an important input to hydrological models. In current hydrological modelling practices, soil property information is often derived from soil category maps by the linking method in which a representative soil property value is linked to each soil polygon. Limited by the area‐class nature of soil category maps, the derived soil property variation is discontinuous and less detailed than high resolution digital terrain or remote sensing data. This research proposed dmSoil, a data‐mining‐based approach to derive continuous and spatially detailed soil property information from soil category maps. First, the soil–environment relationships are extracted through data mining of a soil map. The similarity of the soil at each location to different soil types in the soil map is then estimated using the mined relationships. Prediction of soil property values at each location is made by combining the similarities of the soil at that location to different soil types and the representative soil property values of these soil types. The new approach was applied in the Raffelson Watershed and Pleasant Valley in the Driftless Area of Wisconsin, United States to map soil A horizon texture (in both areas) and depth to soil C horizon (in Pleasant Valley). The property maps from the dmSoil approach capture the spatial gradation and details of soil properties better than those from the linking method. The new approach also shows consistent accuracy improvement at validation points. In addition to the improved performances, the inputs for the dmSoil approach are easy to prepare, and the approach itself is simple to deploy. It provides an effective way to derive better soil property information from soil category maps for hydrological modelling. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Glacial conditions that contribute to the regeneration of Fountain Glacier proglacial icing,Bylot Island,Canada          下载免费PDF全文

Pablo Wainstein  Brian Moorman  Ken Whitehead 《水文研究》2014,28(5):2749-2760

Proglacial icings are one of the most common forms of extrusive ice found in the Canadian Arctic. However, the icing adjacent to Fountain Glacier, Bylot Island, is unique due to its annual cycle of growth and decay, and perennial existence without involving freezing point depression of water due to chemical characteristics. Its regeneration depends on the availability of subglacial water and on the balance between ice accretion and hydro‐thermal erosion. The storage and conduction of the glacial meltwater involved in the accretion of the icing were analyzed by conducting topographic and ground penetrating radar surveys in addition to the modelling of the subglacial drainage network and the thermal characteristics of the glacier base. The reflection power analysis of the geophysical data shows that some areas of the lower ablation zone have a high accumulation of liquid water, particularly beneath the centre part of the glacier along the main supraglacial stream. A dielectric permittivity model of the glacier – sediment interface suggests that a considerable portion of the glacier is warm based; allowing water to flow through unfrozen subglacial sediments towards the proglacial outwash plain. All these glacier‐related characteristics contribute to the annual regeneration of the proglacial icing and allow for portions of the icing to be perennial. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications     

Antonio Sánchez‐Navas  Antonio García‐Casco  Agustín Martín‐Algarra 《地学学报》2014,26(6):477-486

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.  相似文献   

Early solar system aqueous activity: K isotope evidence from Allende     

Yun Jiang  Piers Koefoed  Olga Pravdivtseva  Heng Chen  Chun‐Hui Li  Fang Huang  Li‐Ping Qin  Jia Liu  Kun Wang 《Meteoritics & planetary science》2021,56(1):61-76

The alkali element K is moderately volatile and fluid mobile; thus, it can be influenced by both primary processes (evaporation and recondensation) in the solar nebula and secondary processes (thermal and aqueous alteration) in the parent body. Since these primary and secondary processes would induce different isotopic fractionations, K isotopes could become a potential tracer to distinguish them. Using recently developed methods with improved precision (0.05‰, 95% confidence interval), we systematically measured the K isotopic compositions and major/trace elemental compositions of chondritic components (18 chondrules, 3 CAIs, 2 matrices, and 5 bulks) in the carbonaceous chondrite fall Allende. Among all the components analyzed in this study, CAIs, which formed initially under high‐temperature conditions in the solar nebula and were dominated by nominally K‐free refractory minerals, have the highest K₂O content (average 0.53 wt%) and have K isotope compositions most enriched in heavy isotopes (δ⁴¹K: ?0.30 to ?0.25‰). Such an observation is consistent with previous petrologic studies that show CAIs in Allende have undergone alkali enrichment during metasomatism. In contrast, chondrules contain lower K₂O content (0.003–0.17 wt%) and generally lighter K isotope compositions (δ⁴¹K: ?0.87‰ to ?0.24‰). The matrix and bulks are nearly identical in K₂O content and K isotope compositions (0.02–0.05 wt%; δ⁴¹K: ?0.62 to ? 0.46‰), which are, as expected, right in the middle of CAIs and chondrules. This strongly indicates that most of the chondritic components of Allende suffered aqueous alteration and their K isotopic compositions are the ramification of Allende parent‐body processing instead of primary nebular signatures. Nevertheless, we propose the small K isotope fractionations observed (< 1‰) among Allende components are likely similar to the overall range of K isotopic fractionation that occurred in nebular environment. Furthermore, the K isotope compositions seen in the components of Allende in this study are consistent with MC‐ICP‐MS analyses of the components in ordinary chondrites, which also show an absence of large (10‰) isotope fractionations. This is not expected as evaporation experiments in nebular conditions suggest there should be large K isotopic fractionations. Nevertheless, possible nebular processes such as chondrules back exchanging with ambient gas when they formed could explain this lack of large K isotopic variation.  相似文献